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Abstract 


The vibrational fundamental modes of 7-methoxy- 4-bromomethylcoumarin (7BMC) have been analyzed 
by combining FTIR, FT-Raman and quantum chemical calculations. The structural parameters of the compound 
are determined from the optimized geometry by B3LYP method with 6-311++G (d, p) basis set. The harmonic 
vibrational frequencies were calculated and the scaled values have been compared with experimental FTIR and 
FT-Raman spectra. 'H and °C NMR spectra have been analyzed and 'H and "°C nuclear magnetic resonance 
chemical shifts are calculated using the gauge independent atomic orbital (GIAO) method. The theoretical 
UV—-VIS spectrum of the compound and the electronic properties, such as HOMO (highest occupied molecular 
orbital) and LUMO (lowest occupied molecular orbital) energies were performed by time-dependent density 
functional theory (TD-DFT) approach. Information about the size, shape, charge density distribution and site 
of chemical reactivity of the molecule have been obtained by mapping electron density isosurface with molecular 
electrostatic potential (MESP). The change in electron density (ED) in the o* antibonding orbitals and 
stabilization energies E(2) have been calculated by natural bond (NBO) analysis to give clear evidence of 


stabilization originating in the hyper conjugation of hydrogen—bonded interactions. 
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1. INTRODUCTION 


Coumarins are an important class of oxygen 
heterocycles, widespread in the plant kingdom (Borges 
et al. 2005). Coumarins are of great interest due to their 
biological properties. In particular, their physiological, 
bacteriostatic and anti-tumour activities make these 
compounds more attractive in the field of medicine. 
They have been shown to be useful as antitumoural 
and anti-HIV agents. Some coumarins emit a strong 
fluorescence, on the substitution of various functional 
groups at different positions, with a pronounced charge 
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transfer (CT) character (Prasad et al. 2013). 
4-methyl coumarins have been found to possess 
choleretic, analgesic, anti-spermatogenic, anti- 
tubercular and diuretic properties. 7-methoxycoumarin 
naturally occurs in chamomile and presents a alsamic 
sweet odor. 4-methyl-7-methoxycoumarin is used as a 
fluorescent label for high performance liquid 
chromatography (Prabavathi and Senthil Nayaki, 2014). 
4-bromomethyl-7-methoxycoumarin is used as a 
fluorescent label for fatty acids. It has been used in the 
flurometric analysis of prostaglands and a wide range 
of naturally occurring acids, including bile and 
thromboxane B2 (Alekseev, 1981). Owing to the 
importance listed above, a comprehensive DFT study 
has been performed on the titled molecule. 
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2. EXPERIMENTAL DETAILS 


The compound under investigation 
7-methoxy-4-bromomethyl coumarin (7BMC) was 
purchased from the Sigma Aldrich Company and used 
as such without any further purification. The room 
temperature FT-IR spectrum of 7BMC was measured in 
the region 4000-400 cm"! on BRUKER IFS 66V 
spectrometer usingKBr pellet with the spectral 
resolution +1 cm". The FT-Raman spectrum of the 
chosen molecule was recorded in the region 
4000-100 cm! with BRUKER RFS 100/s Raman module, 
equipped with an Nd:YAG laser source, operating at 
1064 nm linewidth 150 mW power. The reported 
wavenumbers are expected to be accurate to within 
+] cm". 


'H and °C nuclear magnetic resonance (NMR) 
(400 MHz; CDCI,) spectra were recorded on a Bruker 
HC400 instrument. Chemical shifts are reported 1n parts 
per million scale (delta scale) downfield from 
tetramethylsilane. 


3. COMPUTATIONAL DETAILS 


The quantum chemical calculation of 7 BMC 
has been performed using Gaussian 09 program (Benito 
Reyes-Trejo et al. 2014) at the B3LYP level, 
supplemented with the standard 6-3 11++G** basis set. 
The optimized geometries corresponding to the 
minimum potential energy surface have been obtained 
by solving self-consistent field equation iteratively. 
Harmonic vibrational wave numbers have been 
calculated using analytic second derivatives to confirm 
the convergence to minima on the potential energy 
surface and to evaluate the zero-point vibrational 
energy. Multiple scaling of the force field has been 
performed by the scaled quantum mechanical (SQM) 
procedure, to offset the systematic errors caused by 
basis set incompleteness, neglect of electron correlation 
and vibrational anharmonicity. Normal coordinate 
analysis on 7BMC has been performed to obtain full 
description of the molecular motion pertaining to the 
normal modes using the MOLVIB-7.0 program. 


The calculated analytic force constants are used by 
MOLVIB in the calculation of vibrational frequencies 
by diagonalization of dynamical matrix. These force 
fields obtained in cartesian coordinates and dipole 
derivatives with respect to atomic displacements were 
extracted from the archive section of the Gaussian 09 
output and are transformed to a suitably defined set of 
internal coordinates (the ‘natural coordinates’). The 
Raman activities (S;) calculated by Gaussian 09 program 
have been suitably adjusted by the scaling procedure 
of MOLVIB and subsequently converted to relative 
Raman intensities (J) using the following relationship 
derived from the basic theory of Raman scattering. 


4 
ff Wy-V;) S; 
i v,[1-exp(—hev, / kT’) 


where v, is the exciting frequencuy(in cm” units), v, is 
the vibrational wave number of ther ith normal mode, h, 
c, k are the fundamental constants and f is a suitably 
chosen common normalization factor for all peak 
intensities. 


The symmetry of the molecule is also helpful 
in making vibrational assignments. By combining the 
results of the Gauss view program with symmetry 
considerations, vibrational frequency assignments are 
made with a high degree of confidence. There is always 
some ambiguity in defining internal coordinates. 
However, the defined coordinates form a complete set 
and matches quite well with the motions observed using 
the Gauss view program. 


The natural bonding orbital (NBO) calculations 
were performed using NBO 3.1 program, implemented 
in the Gaussian 09 package at the DFT/B3LYP level. 
NBO helps to understand various second order 
interactions between the filled orbitals of one subsystem 
and the vacant orbitals of another subsystem, which is 
a measure of the intermolecular delocalization and hyper 
conjugation. The hyper conjugative interaction energy 
was deduced from the second-order perturbation 
approach (Prabavathi et al. 2013). 
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(2) F? 
E = AE = q 


where Fi; is the Fock matrix element between 1 and j 
NBO orbitals ¢,and € are the energies of o and o* NBO’s, 
and q, is the population of the donor o orbital. 


4. RESULTS & DISCUSSION 
4.1 Molecular Geometry 


In order to figure out the minimum energy 
conformer of the titled molecule, potential energy 
surface (PES) scan has been performed on the dihedral 
angle Br16-C13-C4-C3. Structure optimization is done 
for the minimum energy conformer. The optimized 
structure of 7BMC with numbering scheme for the atoms 
is presented in Fig.1 and it belongs to C, point group 
symmetry with energy of about -3224.583 Hartrees. The 
optimized structural parameters like bond length, bond 
angle and the dihedral angle of 7BMC determined by 
B3LYP method with 6—311++G** basis sets are 
summarized in the Table 1. All the C-H bond length of 
TBMC averages to 1.08 A except methoxy and 
bromomethyl groups, where it averages to 1.09 A. In 
TBMC the O11-C2bond length is 1.2 A and it coincides 
with the reported C-O bond length 1.2 A of 7-methoxy- 
4-methyl coumarin, 3-acetylcoumarin and 3-cyano 
methyl coumarin (Prabavathi and Senthil Nayaki, 2014). 
The ring C-O bond length is usually 1.43 A (Arjunan 
et al. 2013). But, a reduction in bond length (C2-Ol= 
1.39 and C9-O1=1.36 A) observed in the present 
investigation is due to the fusion of benzene ring with 
the a—pyrone ring through O1 atom (Prabavathi and 
Senthil Nayaki, 2014). The substituent plays a vital role 
on the structural and electronic properties of the 
molecules. The bromomethyl and methoxy group 
substitution affect the regular hexagonal structure of 
their respective rings, which is evident from the 
observed bond angles. The decrease in bond angle C7- 
C8-C9 (118°) in turn reduces the C7-O19 bond length 
(1.35 A). 


4.2 Hyperpolarizability 


The first hyperpolarizability (Bq) and its 
components of the 7BMC, calculated based on the 


finite-field approach are given in Table 2. The ground- 
state dipole moment and the hyperpolarizability are 
governed by the extent of z electron delocalization, 
which in turn depends on the structural details of the 
molecules. In the current work the dipole moments, as 
well as first hyperpolarizabilities are nonzero as the 
equilibrium geometries are noncentrosymmetric. The 
dominant direction of charge delocalization is found to 
be along B,. Asthe B,, value of 7BMC is greater 
than urea (Sebastain et al. 2011), this can serve as a 
good nonlinear optical material [NLO]. 


4.3 Analysis of Molecular Electrostatic Surface 
Potential (MESP) 


Molecular electrostatic potential (MESP) 
mapping is very useful in the investigation of the 
molecular structure with its physiochemical property 
relationships (Arjunan et al. 2014). The MESP has 
been also employed as an information tool of chemistry 
to describe different physical and chemical features, 
including non-covalent interactions in complex 
biological system. It provides a visual method to 
understand the relative polarity of the molecule. The 
different values of the electrostatic potential are 
represented by different colors; red represents the 
negative regions of electrostatic potential, blue 
represents the positive regions of electrostatic potential 





Fig. 1: Molecular structure of 7-methoxy-4- 
bromomethylcoumarin with atom numbering scheme. 
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Table 1. Optimized geometrical parameters of 7-methoxy-4-bromomethylcoumarin. 


Bond length Value (A ) Bond angle 
C2-01 1.39 C3-C2-O1 
C3-C2 1.45 C4-C3-C2 
C4-C3 1.36 C10-C4-C3 
C10-C4 1.45 C5-C10-C4 
C5-C10 1.41 C6-C5-C10 
C6-C5 1.38 C7-C6-C5 
C7-C6 1.41 C8-C7-C6 
C8-C7 1.39 C9-C8-C7 
C9-C8 1.39 O11-C2-O1 
O11-C2 1.20 H12-C3-C2 
H12-C3 1.08 C13-C4-C3 
C13-C4 1.50 H14-C13-C4 

H14-C13 1.09 H15-C13-C4 
H15-C13 1.09 Br16-C13-C4 
Brl16-C13 1.98 H17-C5-C10 
H17-C5 1.08 H18-C6-C5 
H18-C6 1.08 019-C7-C6 
O19-C7 1.35 C20-019-C7 
C20-O19 1.43 H21-C20-019 
H21-C20 1.09 H22-C20-0 19 
H22-C20 1.09 H23-C20-0 19 
H23-C20 1.09 H24-C8-C7 
H24-C8 1.08 











Value ( °) Torsional angle Value (°c ) 
115.72 C4-C3-C2-O1 0 
122.97 C10-C4-C3-C2 0 
119.28 C5-C10-C4-C3 180 
125.52 C6-C5-C10-C4 180 
121355 C7-C6-C5-C10 0 
120.01 C8 -C7-C6-C5 0 
120.01 C9-C8-C7-C6 0 
118.93 O11-C2-O1-C9 180 
117.92 H12-C3-C2-O1 180 
115.02 C13-C4-C3-C2 180 
118.01 H14-C13-C4-C3 0 
110.57 H15-C13-C4-C3 -120 
111.70 Br16-C13-C4-C3 120 
114.97 H17-C5-C10-C4 0 
119.6 H18-C6-C5-C10 180 
121.31 O19-C7-C6-C5 180 
113.7 C20-019-C7-C6 180 
118.89 H21-C20-019-C7 180 
105.69 H22 -C20-0O19-C7 61 
111.06 H23 -C20-0O19-C7 -61 
111.27 H24-C9-C8-C6 180 
122.98 














Fig. 2: The total electron density isosurface mapped 
with molecular electrostatic potential of 7-methoxy- 
4-bromomethyl coumarin 


and green is the region of zero potential. Potential 
increases in the order red < orange < yellow < green < 
blue (Sheela et al. 2014). The projection of such a 
mapped electrostatic potential surface for the title 
molecule is given in Fig. 2. In the present work, the 
calculated result shows that the negative potentials are 
mainly over the electro negative oxygen atoms O1 and 


O11. The dominance of blue colour over the hydrogen 
atoms of 7BMC indicates the presence of positive 
potential near hydrogen atoms. Thus the hydrogen 
atoms are the preferred sites for nucleophilic reactivity. 
The extensive delocalization of charges from the ring 
carbon atoms and bromine atom is proved by the 
prevalence of green colour over these atoms. 


4.4 Thermodynamic Properties 


The zero point energies, thermal correction to 
internal energy, enthalpy, Gibbs free energy and entropy 
and heat capacity for a molecular system were computed 
from the Gaussian frequency calculations. The computed 
thermodynamic parameters are listed in Table 3. On the 
basis of vibrational analysis, the thermodynamic 
functions such as heat capacity, entropy (S) and 
enthalpy change (AH, ,,) at constant pressure (Cp), for 
7BMC were also determined and are listed in Table 4. 
The graph showing the correlation of heat capacity, 
entropy (S) and enthalpy change (AH, ,,) with 
temperature at constant pressure (Cp), 1s delineated in 
Fig. 3. All the computed thermodynamic functions of 
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Table 2. Calculated all £ components and £6 


TOT 


(x10-30) value of 7-methoxy-bromomethylcoumarin 











B3L Y P/6-311++G** 
B -components 7-methoxy-4- 
bromomethylcoumarin 

Deg -804.52 
Die 521.46 
— -334.93 
Byyy -124.50 
Baxz 10.25 
Byyz -46.85 
Byyz 43.20 

— -6.80 
By oe -32.49 
Dea 6.71 

Beotai(€su) 10.4*10°° 
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Fig. 3: Temperature dependence of heat capacity, 
Entropy and Enthalpy change at constant pressure of 
7-methoxy-4-bromomethylcoumarin 


7BMC increases steadily with temperature, ranging from 
100 to 1000 K, due to the enhancement of the molecular 
vibration with increase in temperature. 


Here, all the mentioned thermodynamic 
calculations were done in gas phase. As per the second 
law of thermodynamics in thermo chemical field 
(Sebastian et al. 2013), these calculations can be used 
to compute the other thermodynamic energies and helps 
to estimate the directions of chemical reactions. 


4.5 Spectral Analysis 


Detailed description of vibrational modes can 
be given by means of normal coordinate analysis (NCA). 
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Fig. 4: Comparison of observed and calculated FT- 
RAMAN spectra of 7-methoxy-4-bromomethyl 
coumarin (a) observed and (b) calculated with B3LYP/ 
6-311++G** 
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Fig. 5: Comparison of observed and calculated FTIR 
spectra of 7-methoxy-4-bromomethylcoumarin (a) 
observed and (b) calculated with B3LY P/6-311++G** 


For this purpose, the standard internal coordinates 
(containing redundancies) were constructed. From 
these, nonredundant set of local symmetry coordinates 
were constructed by suitable linear combinations of 
internal coordinates Elbert et al. (1985). In agreement 
with the Cl symmetry, the 66 normal modes of 7BMC 
are distributed among the symmetry species as 


T3n_6= 45 (in-plane) + 21 (out-of-plane) 


The detailed vibrational assignment of 
fundamental modes of the titled molecule was made 
unambiguously with the help of Gauss view program 
and the PED (potential energy distribution) obtained 
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Table 3. The calculated thermodynamic parameters of 7-methoxy-4-bromomethylcoumarin. 


Structural parameters 


SCF energy (Hartrees) 


Total energy (thermal), Etota] (kcal mol!) 


Heat capacity at const. volume, Cv (cal mol-!K-!) 
Entropy, S (cal mol-lK-!) 
Vibrational energy , Eyjh (kcalmol 1) 
Zero-point vibrational energy, Eg (kcal mol!) 
Rotational constants (GHz) 


A 
B 
C 


Dipole moment (Debye) 


Ux 

By 

Uz 
Mtotal 


SCF GIAO Magnetic shielding 


Degeneracy 











80 
Shift (ppm) 


6-311++G** 
-3224.5831 
119.658 
49.573 
117.043 
117.881 


111.382 


0.58843 
0.28876 
0.20058 





SCF GIAO Magnetic shielding 


l4H 23H 





18H 24H = = =12H 


Degeneracy 








$a 40 


Shift (ppm) 





Fig. 6: Theoretically simulated NMR spectrum of 13C and 1H of 7-methoxy-4-bromomethyl coumarin 


during the scaling procedure with Molvib Program. 
These assignments along with calculated IR, Raman 
intensities and normal mode description (characterized 
by PED) are reported in Table 5. The observed and 
simulated FI-Raman and FT-IR spectra of 7BMC were 
produced in common frequency scales and are 
presented in Figs. 4 and 5 respectively. The RMS error 
between unscaled and experimental frequencies was 
reduced to 4.7 by applying refinement of scaling factors 
and optimizing via least square refinement algorithm. 


4.5.1 C-C vibrations 


In the present work, the experimental and 
theoretical skeletal (ring CC) vibrations were assigned 
in the region 915-1630 cm" (Prabavathi and 
Senthilnayaki, 2014). The peaks observed at 1621, 1607, 
1557, 1400, 1357, 1293, 1214, 1057 and 917 cm! in FI-IR 
spectrum are assigned to CC-stretching vibrations. 
Except three all these bands were present in FI-Raman 
spectrum. The experimental data shows excellent 
agreement with the theoretical data. The CC in-plane 
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Fig. 7: The linear regression between the experimental and theoretical 'H and “C NMR chemical shifts of 7- 
methoxy-4-bromomethylcoumarin. 
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Fig. 8: Illustration of Natural atomic charges of 7- 


methoxy-4-bromomethylcoumarin Fig. 9: Natural atomic charges analysis chart 


Table 4. The temperature dependence of 
thermodynamic parameters of 7-methoxy-4- 
bromomethylcoumarin 


LUMO PLOT 








Entro Heat capacit Enthal EACLE 
Temp. (K) al mol K) ‘(Om Pe K} (KJ not) Wes 2 ean veieveN 
100 33145 96.14 6.27 2S, ”* 
200 41640 155.90 18.86 AE = 2.9452 ev 
300 491.16 216.84 37.51 re $. 
400 561.47 273.47 62.08 Epouc™ -8-9669 eV 
500 | 62781 321.27 91.90 a. ee 
600 689.93 259.9] 126.03 eens 
700 74785 391.17 163.65 -4a> ~ 
800 801.80 416.60 204.08 
000 95012 437.60 246 82 se af The shone ya ee a of 
1000 899.16 455.16 29] 48 rontier molecular orbital for /-methoxy-4- 





bromomethylcoumarin 
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Table 5. Observed and B3LYP / 6-311++G** level calculated vibrational frequencies (in cm “D of 7-methoxy-4- 


OMmAAINDNNBWN 


Observed frequencies 


Gree 


Infrared Raman 


Calculated frequencies (cm 


Unscaled 


bromomethylcoumarin 
-1) with B3LYP/ 


6-311++G** level 


Scaled 





TRaA,) 


61.31 
1.521 
90.06 
41.68 
22169 
523.01 
646.00 
528.54 
40.626 
703.098 
584.229 
261.143 
141.656 
166.592 
3.127 
11.876 
176.54 
100.132 
520.05 
118.94 
197.649 
474.76 
397.87 
2.331 
60.907 
182.75 
7.873 
11.242 
196.20 
398.12 
493.49 
286.740 
2.183 
144.57 
S122 
66. 84 
PIES 
245.76 
343.60 
214.45 
112.807 
91.16 
311.01 
0.214 
11.68 
21.31 
104.27 
204.44 
125.26 
63.68 
122.27 
43.70 
ZlA23 
2.801 
0.853 


Raman 


76.400 
64.391 
111.830 
117.562 
50.766 
172.524 
97.144 
51.425 
149.173 
173.554 
105.275 
48.412 
254.462 
81.078 
20.754 
6.225 
15.690 
92121 
16.042 
84.338 
78.072 
65.625 
4.000 
102.620 
97.010 
86.517 
98.261 
4.125 
3.703 
0.321 
34.428 
23.974 
17.149 
17.678 
0.288 
1.837 
3.641 
5997 
0.254 
0.875 
17.812 
3.993 
6.690 
7.656 
7.419 
8.670 
2.293 
6.454 
23.019 
10.598 
4.970 
5.086 
4.891 
21.905 
2.701 


bil.) 
i 


PED (% ) among type of 
internal coordinates 


CH(99) 

CH(97) 

CH(99) 

CH(97) 

CH2ASS(80),CH2SS(19) 

CH2SS (62),CH2ASS(19), 

OCH3o0p(79), OCH3ip (20) 
OCH3ip(78),OCH3o0p(21) 
OCH3SS(80),O0C H30p(17) 

CO(80),brin g1(9),CC(7) 
CC(32),C03(18),bCCH(15),CC (14) 
CC(29),CC(22),CC(18), bring2(6) 
CC(34),CC(27) bCCH(10),CC25) 
CC(32),bCCH(22),CC(11),CO(10) 

OCH 3ipb(88),OCH3ips(5) 

OCH 3opb(81),O0CH3ops(13) 
CH2SCI(94) 

OCH3sb(70), bCCH(10) 
CC(37),bCCH(17),CC2(15) 
bCCH(18),CC(14), bring] (10),CC (10) 

CC (25),CO03(20),bCCH(9),OCH3sb(15) 
CC(30),CC3(17),CC1 (13),bCCH(12) 
bCCH(65),CC(8),CH2TW(6),CC(6) 
CH2TW (41),CC(24),bring1(12), bCCH (9) 
bCCH(22),CH2ROC(13),CC(13), bring2(10) 
CC (23),bCCH(16),CO2(13) 
CH2ROC(37),OCH3 o0p(23),bCCH(16),CC4(6) 
OCH 3ipr (22),CO2(15),CH2ROC(15),bCCH(12) 
OCH 3ipr(80),OC H30pr( 14) 
bCCH(51),CC(21),CO(Q), bring] (7) 
CC(25),CO (19),CC2(13),bring1(11) 
CO(35),CO(16),CC(12),CC(11) 
CO(23),CCU9),CC(18),C C14) 
CO(39),CO(20), bring2(8),CC(6) 

gCH(8 9), tring2(10) 

CC(51),gCH(12), CH2ROC(Q) 
gCH(68),gCO(9),tring1(7) 
CO(34),tCH2(13),CC(9),gCH (6) 
g¢CH(66),tring2(23),gO0CO(7) 
gCH(72),tring2(10),tring 1(6) 
bring1(36),CC2(18), bring2(14) 
bring1(38),bring2(31),g¢CO(10) 
bOC(32), tring 1(7),CBr(7),gCC(7) 
bOC(19),bring2(15),tring2(14) ,CC(9) 
bring 1 (33),bring2(14),CC2(8) 
CBr(17),bring2(14),bCCH(14) ,bCBr(10) 
bCOC(38),tring2(34),tring1(11), gC H(6) 
bCBr(30),bring2(19),bOC(19),CC1(8) 
bring2(25),tring! (14), bring 1(10) 
bCC(15),tButt(14),tring1(11),bring1(10) 
gCOC(36),bring1(21),bring2(10) ,bOC(8),bCO(7) 
bring2(45), g0CO(13),tButt(10) 

gOC(37), bOC(5)bring2(12),tring1(5), 
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Observed Calculated frequencies (cm 


S.No. 


Infrared Raman Unscaled Scaled IR@ 


57 
58 
a9 
60 
61 
62 


63 
64 








-1) with B3LYP/ 


frequencies(em~)) 6-311++G ** level 


(Aj) 


1,359 
2.507 
1.048 
2212 
2.841 
1.184 


0.121 


L925 
2.124 
1.486 


Abbreviations used: b, bending; g, deformation; t, torsion; Butt, butterfly. 


a Relative absorption intensities normalized with highest peak 


absorption. 


b Relative Raman intensities calculated and normalized to 100. 


and out-of-plane bending vibrations les in the region 
1000-400 cm! (Prasad et al. 2013). So, the band observed 
at 507 cm!‘ in FT-IR and at 508 cm" in FI-Raman spectra 
is readily assigned for CC in-plane bending vibration. 
The out-of- plane vibration 1s allotted for the band at 
340 cm" in FI-Raman spectrum. 


4.5.2 C-H vibrations 


The characteristic C—H stretching vibrations 
of heteroaromatic structure are expected to appear in 
3100-3000 cm" frequency range. In infrared spectra, 
most of the aromatic compounds have nearly four peaks 
in the region 3080-3010 cm" due to ring C—H stretching 
bands. IR frequencies of C—H bands are a function of 
sp hybridization (Prabavathi et al. 2014). Among the 
four C-H stretching vibrations, three modes are 
identified in IR spectrum at 3010, 3030 and 3085 cm" 
and the fourth in Raman at 3075 cm”. These vibrations 
are in good agreement with the calculated values, with 
an average of 98% PED contribution. 


The C-—H in-plane-bending vibrations are 
usually expected to occur in the region 1300-1000 cm" 
and these vibrations are very useful for characterization 
purposes (Pavia et al. 2001). The calculated values 
1140, 1263, 1271 and 1375 cm‘ are assigned to four C-H 
in-plane-bending vibrations. All the four peaks are 
present in FTIR spectrum at 1129, 1236, 1271 and 
1375 cm’. The C-H out-of-plane bending vibrations 
are strongly coupled and occur in the region 


RamanP (Ii) 





PED (% ) among type of 
internal coordinates 







tring 1(22),tOC H3(19),tCH 2(1 1) 
tring 1(35),tring2(30),gCH(10) 
tring] (23),bCC(22),tCH2(17) 
tOCH3(42),tring 1(29),tring2(10),g0CO(6) 
tring2(24),tCBr(22),tCH2(16) 
tring2(25),bCO(23),tCH2(15) 
tring2(21),tCOC(19),tring 1(18) 


tButt (72),tring2 (8),gCH(8) 


tCOC(44),tOC H3(14), gCC(14) 
tCBr(51),tCH2(42) 
tCH2(20),tCOC(20), eC C(13) 

















Table 6. Experimental and Theoretical isotropic 
chemical shifts (with respect to TMS, all values in 
ppm) for 7-methoxy-4-bromomethylcoumarin 





Chem ical Shift (ppm ) 
Calculated Ex perim ental 
164 156 
160 150 
132 126 
118 114 
116 113 
115 111 
102 101 
56 56 
46 pie | 
8 .1 7.6 
7.1 Teo 
6 .8 6.8 
6.1 6.4 
4.7 4.4 
4.2 3.8 
4.1 4.4 
3.8 3.8 
3.8 3.8 








Table 7. The Second-order perturbation energies E(2) 
(Kcal/mol) corresponding to the most important 


charge transfer interactions (donor-acceptor) in 7- 
methoxy-4-bromomethylcoumarin by 


Excited 


state 


S1 
S2 
S3 
S4 





B3LYP/6-311++G** method 





Excitation Oscillator 
Wavelength : 
Wwerey energies strengths 
(eV) (f) 
338.6 3.66 0.29 
214.3 5.78 0.13 
201.4 6.15 0.42 
200.7 6.17 0.22 





Table 8. The Second-order perturbation energies 
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E(2) (Kcal/mol) corresponding to the most important charge 


transfer interactions (donor-acceptor) in 7-methoxy-4-bromomethylcoumarin by B3LYP/6-311++G** method. 






























UV -VIS Spectrum 
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Fig. 11. Theoretically calculated UV-Vis spectrum of 7- 


methoxy-4-bromomethylcoumarin 











. E(2)4 E(Qj)-E@ F(i,j)©& ED,s(%) ED 3 (% ) 
NBO ED/ NBO ED/ 
@ E (i) 5 (kJmol) b (au.) (au) (aujofj (au) of j 
0.02303 
(c) 0.03493 
01-C2 1.99011 0.01250 
0.02723 
0.01228 
0.01943 
is) 0.02303 
01-C9 1.98957 0.03493 
0.01250 
0.02723 
a 0.01943 
1.98139 0.02303 
Cer? 0.03493 
(o) 
ee 1.99444 0.37600 
C2-O11 0.01228 
C2-011 0.28458 
Pee 1.97452 C3-C4 0.01943 
C20-H21 0.00847 
C20-H22 0.01736 . 
(x) naan C3-C4 0.01943 27.89 0.41 0.091 50.96 49.04 
C3-C4 C10-C9 0.03493 20.01 0.28 0.069 49.15 50.85 
(o) 
ae 1.97516 C2-011 0.28458 14.90 
(o) 
er 1.96776 C10-C5 0.02303 27.54 
(s) Tees C10-C5 0.02303 14.89 0.91 0.104 48.26 51.74 
C10-CS5 C10-C5 0.02303 32.03 0.93 0.155 48.26 51.74 
is) 0.03493 
C10-C9 aad 0.02723 
(c) 0.03493 
C5-C6 ev Ore 0.02723 j 
my 1.96288 C7-C8 0.02723 41.15 0.87 0.169 50.06 49.94 


1000 to 700 cm" (Prabavathi et al. 2014). Accordingly, 
the calculated values are assigned as shown in Table 5.. 


4.5.3 Bromomethyl vibrations 


The C-Br stretching vibration computed at 
667 cm” exactly coincide the experimental value, in both 
IR and Raman spectrum. The vibrational frequencies and 
IR intensities of CH,-Br group are sensitive to 
conformational changes and nature of substituent 
(El-Nahass et al. 2011). The CH,-Br asymmetric 
stretching vibration is identified at 2975 cm” in Raman 
and the findings given in the literature 
(Prabhavathi et al. 2014) supports this assignment. The 
CH,-Br symmetric stretching vibration appears in 
FT-IR and FT-Raman spectra at 2950 and 2942 cm", 
respectively. The fundamental CH,-Br vibrations due to 
scissoring, wagging, twisting and rocking appear in the 
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frequency region 1500-800 cm (Silverstein et al. 1981). 
In light of the above statement the computed values at 
1197, 1269 and 1450 cm" were assigned for CH,-Br 
rocking, twisting and scissoring vibrations, respectively. 
Among the vibrations mentioned above, the rocking 
mode vibration is present in IR spectrum at 1193 cm" 
and the other two vibrations were found in Raman at 
1266 and 1450 cm" respectively. 


4.5.4 C-O Vibrations 


The C=O group vibrations are likely to be the 
most sensitive to the environment. Hence, they show 
pronounced shifts in the spectra of the aromatic 
compounds. The C=O stretching band is characterized 
by a very strong and sharp band appearing 
approximately at 1700 cm'(Balachandran and 
Karunakaran, 2013). So, for 7BMC the C2=O11 
stretching vibration is allotted for the simulated value 
1720 cm and its experimental counterpart in IR spectrum 
exactly coincide with this value.The two ring CO 
stretching vibrations (C9-O1 and C2-O1) are predicted 
at 851 and 983 cm''. Both of these vibrations are found 
in IR at 865 and 990 cm", respectively. The methoxy 
group C-O stretching vibrations are calculated at 1024 
and 1037 cm'1. Out of these twobands, one is found in 
IR spectrum alone at 1029 cm” and the other is present 
in Raman spectrum only at 1025 cm”. 


4.5.5 Methoxy group vibrations 


The compound selected for the present study 
has one methoxy group substituted in the 7" position 
of the titled molecule. For the assignments of CH, 
group frequencies, one can expect that nine 
fundamentals can be associated to CH, group namely 
CH,ss, symmetric stretch; CH, ips, in-plane stretch (i.e. 
in-plane hydrogen stretching mode); CH, ipb, in-plane 
bending, (i.e. hydrogen deformation mode); CH,sb, 
symmetric bending; CH, ipr, in-plane rocking; CH, opr, 
out-of plane rocking; CH, ops, out-of-plane stretch; 
CH, opb, out-of plane bending modes and tCH, twisting 
modes. Among which the CH, ops, out-of-plane stretch 
and CH, opb, out-of plane bending modes of CH, group 
would be expected to be depolarized for A” symmetry 
species. In addition to this OC wagging and C-O—C 


in-plane bending also exist. The CH, symmetric and 
asymmetric stretching vibrations generally occur at 
lower frequency region 2850-3000 cm'!. The out-of- 
plane stretching mode of CH3 group is expected around 
2980 cm! and symmetric one is expected around 
2870 cm! (Krishnakumar and Prabhavathi, 2010). For 
the title compound the peaks identified in IR at 2920 
and 2860 cm" and in FI-Raman at 2850 cm‘ are attributed 
to CH, ops, CH, ips, and CH, ss vibrations respectively. 


For methyl substituted benzene derivatives the 
antisymmetric and symmetric deformation vibrations of 
methyl group normally appear in the region 
1465-1440 cm" and 1390-1370 cm" respectively and 
the rocking vibrations of CH3 modes usually appear in 
the region 1070-1010cm" (Krishnakumar and 
Prabhavathi, 2010). Based on the above data the bands 
found in FIT-Raman at 1150,1183,1433, 1458 and 
1517 cm" are assigned for the vibrations OCH,opr, 
OCH, ipr, OCH,sb, OCH, opb and OCH, ipb, respectively. 
The torsional vibration (tOCH3) evaluated at 221 cm! 
was found in Raman spectrum at 216 cm’. The 
experimental values agree well the simulated values. 


4.6 NMR Analysis 


The experimental values for 'H and °C NMR 
of 7BMC were measured in CDCL,solution. The 'H and 
SC. theoretical and experimental chemical shifts are 
tabulated in Table 6. The correlation coefficient for 
proton (0.991) and carbon (0.993) chemical shift is 
calculated from linear regression curves. These values 
confirm the good agreement between experimental and 
calculated chemical shift values. The theoretical '°C 
and 'H NMR spectra of 7BMC and their respective linear 
regression curves were depicted in Fig. 6 and Fig. 7, 
respectively. Hydrogen attached to near by electron 
withdrawing atom or group can decrease the shielding 
and move the resonance of attached proton towards a 
higher frequency. By contrast electron donating atom 
or group increases the shielding and moves the 
resonance towards to a lower frequency (Sebastin 
etal. 2011). Accordingly for 7BMC the chemical shift 
value of —CH,Br and —OCH, attached H atom 
experiences a lesser value of about 3.8 to 4.4 ppm and 
fairly agrees with the calculated value. The other four 
hydrogen atoms attached to ring carbons have a greater 
value. 
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Aromatic carbon gives signals with chemical 
shift values from 100 to 200 ppm (Prabhavathi and Senthil 
Nayaki, 2014). Except C13 and C20 atoms, the chemical 
shift values of all the carbons of 7BMC are in accordance 
with the above statement. The shielding around C13 
atom in bromomethyl group and C20 atom in methoxy 
group of 7BMC is somewhat high enough compared to 
other carbon atoms due to the presence of nearby 
electronegative atoms. So the chemical shifts of those 
carbons are set upfield with a low ppm value of about 
27 and 56 ppm respectively. 


4.7 UV-Vis Spectral Analysis 


The time-dependent density functional theory 
(TD-DFT) calculation has been performed to investigate 
the electronic absorption properties of 7BMC. The 
calculated visible absorption maxima (A__) which are 
the functions of electron availability, theoretical 
electronic excitation energies and oscilla-tor strength 
are all tabulated in Table 7. The theoretically calculated 
UV-VIS spectrum 1s depicted in Fig. 11. The calculated 
Ri values for 7BMC are 338.6, 214.3, 201.4 and 200.7nm, 
among which very strong light absorptions at 214.3, 
201.4 and 200.7 nm are due to x a* transitions. The 
absorption band of 7BMC at the longer wave length 
region 338.6 nm is caused by the n—a* transition. This 
type of transition is the characteristic nature of 
compounds containing double bonds involving hetero 
atoms (e.g., C=O). 


4.8 Natural Atomic Charges 


The calculation of effective atomic charges 
plays an important role in the application of quantum 
mechanical calculation to molecular systems. The 
calculated natural atomic charge values from the natural 
population analysis (NPA) are portrayed in the Fig. 8. 
Natural atomic charges analysis chart is shown in 
Fig. 9. From the figures, it 1s evident that there is no 
difference in charge distribution observed on all the 
hydrogen atoms. The methoxy group hydrogen atoms 
show a slightly less value which is not appreciable. The 
carbon atoms C2, C7 and C9 exhibit positive charges 
due to the presence of electronegative oxygen atoms in 
their neighbourhood. The neutral charge on C4 and Br16 
atom shows the pinnacle delocalization of their charges 
towards C14 atom. So, C14 1s more negative than the 


other carbon atoms of 7BMC. The charge on the oxygen 
atoms of 7BMC averages to -0.5. 


4.9 Frontier Molecular Orbitals 


Frontier molecular orbitals (FMOs) play an 
important role in the optical and electronic properties, 
as well as in quantum chemistry and UV-Vis spectra. 
The HOMO represents the ability to donate an electron 
where LUMO, as electron acceptor, represents the ability 
to obtain electron and the energy gap between HOMO 
and LUMO characterizes the molecular chemical 
stability, chemical reactivity, hardness and softness of 
the molecule (Sylystree et al. 2012). The three 
dimensional plots of the FMOs are shown in Fig. 10. 


The highest occupied molecular orbitals are 
localized mainly on the parent ring and the methoxy 
group. The LUMO plot is almost uniformly distributed 
except the methoxy group. The HOMO and LUMO 
transition implies an electron density transfer from the 
coumarin ring to the bromomethyl group. The transition 
from HOMO to LUMO of the molecule results in the 
energy gap (AE) of about 2.945 eV. From the energy gap 
value it can be shown that the titled compound has low 
band gap energy. Thus the molecule 7BMC implies high 
chemical reactivity, because it is energetically favorable 
to add electron to a high-lying LUMO by extracting 
electrons from low-lying HOMO (Aihara, 1999). 


4.10 NBO Analysis 


Table 8 depicts the bonding concepts such as 
type of bond orbital their occupancies, the natural atomic 
hybrids and the percentage of the NBO on each hybrid 
of 7BMC molecule determined by B3LY P/6-311++G** 
method. The occupancies of NBO’s in 7BMC reflect their 
exquisite dependence on the chemical environment. The 
NBO energy values show the corresponding spatial 
symmetry breaking in the direction of unpaired spin. 
The most important bond—antibonds and lone-pair— 
antibonds interaction energies of 7BMC are given in 
Table 8. The intra-molecular interaction 1s formed by the 
orbital overlap between o (C—C) and o* (C—C) bond 
orbital which results intra-molecular charge transfer (ICT) 
causing stabilization of the system. These interactions 
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are observed as increase in electron density (ED) in 
C-C antibonding orbital that weakens the respective 
bonds.The strongest delocalization in 7BMC is between 
o (C2-O11)— (C7-C8) and o (O1-C9) 4o0*(C10-C9). 
These interactions stabilize the titled molecule by 740.7 
kJ/mol and 126.8 kJ/mol, respectively. The other 
interactions producing a large stabilizing energies are 
6 (O1-C9) o0* (C10-C5), c (O1-C9)— o* (C3-C4) and 
6 (C2-C3) o* (C3-C4), with energies of about 66.84, 
51.22& 43.7 kJ/mol respectively. 


5. CONCLUSION 


In this study, the optimized molecular 
structure, PED, thermodynamicand electronic 
properties of the title compound are calculatedby DFT 
method using B3LYP/6-311G(d,p) basis set. 
Thecomplete molecular structural parameters and 
thermodynamicproperties of the compound have been 
calculated. The correlations between the statistical 
thermodynamics and temperature are also obtained. It 
was seen that the entropies and enthalpies increase 
with the increase in temperature as the intensities are 
enhanced by the higher temperature values. The 
fundamental vibrationalmodes of the title compound 
have been precisely assigned, analyzed and the 
experimental results were compared withthe theoretical 
values. As the B value of the compound is high, 7BMC 
can be used as a good NLO material. The molecule 
7BMC shows high chemical reactivity which is evident 
from the low band gap energy. The influences of 
electronegative atoms on the atomic charges 
werestudied in detail. The prediction of reactive 
behaviorof 7BMC in both electrophilic and nucleophilic 
reactions has been done with the help of MESP 
visualization. The chemical shift shielding tensors were 
predicted by the computational NMR method. Thus 
the complete vibrational assignments, structural 
information and electronic properties of the titled 
compound were provided in the present investigation. 
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